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Abstract 

An exactly scale-invariant spectrum of scalar perturbation generated during de 
Sitter spacetime is found from the gravity model of the nonminimal derivative coupling 
with fourth-order term. The nonminimal derivative coupling term generates a healthy 
(ghost-free) fourth-order derivative term, while the fourth-order term provides an 
unhealthy (ghost) fourth-order derivative term. The Harrison-Zel’dovich spectrum 
obtained from Fourier transforming the fourth-order propagator in de Sitter space 
is recovered by computing the power spectrum in its momentum space directly. It 
shows that this model provides a truly scale-invariant spectrum, in addition to the 
Lee-Wick scalar theory. 
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1 Introduction 


In the evolutionary phase of density inhomogeneities, one often makes the simplifying as¬ 
sumption that the primordial power spectrum has the simple power-law expression like 
V(k) = A s (k / k^ ns ~ l +^r ln (. k l k *) with fc* a pivot scale. The case of n s = 1 and a s = 0 corre¬ 
sponds to the Harrison-Zeld’ovich (HZ) spectrum [I] E3 3] and this has been ruled out by 
different datasets |4j. In the picture of slow-roll inflation (quasi-de Sitter expansion), the 
momentum dependence at any given time arises as a consequence of the time dependence 
of H and 0 compared to the de Sitter (dS) expansion of constant H and 0 = 0. Recent data 
from Planck has shown that the scalar spectrum is a nearly scale-invariant one with the 
amplitude A s = . ^ (H 2 /2tt) 2 = (2.441 ± 0.092) x 10~ 9 , implying that it is approximately 
l(more precisely, n s = 0.9603 ± 0.0073) [5j. 

On the other hand, it is worth noting that the power spectrum of a massless minimally 
coupled (nunc) scalar in dS expansion takes the form of (H/2tt) 2 [1 + ( k/aH ) 2 ]. It reduces 
to the HZ scale-invariant spectrum of (H/2ti) 2 in the superhorizon region of k <C aH , 
whereas it leads to ( k/2na ) 2 in the subhorizon region of > aH [6|. We note here that 
the latter is just the spectrum of a massless conformally coupled (mcc) scalar generated 
during dS expansion. This happens to the tensor spectrum too. Clearly, the HZ scale- 
invariant spectrum is related to the two-point function which is logarithmically growing for 
largely separated points in dS space. Recently, it was shown that the IR growing of the 
two-point function is physical because Fourier transforming logarithmic two-point function 
can lead to the HZ spectrum when one tames logarithmic divergence by using Cesaro- 
summability technique [7J. This implies that one could obtain the HZ spectrum through 
an IR regularization procedure. 

More recently, the authors have obtained the HZ spectrum (H/2n) 2 from the Lee- 
Wick model of a fourth-order derivative scalar theory in dS spacetime [8]. Here we have 
obtained the fourth-order propagator as the inverse of the Lee-Wick (LW) operator Alw = 
— jj 2 (V 4 — M 2 V 2 ) with M 2 mass parameter. The operator of M 2 Alw” becomes the Weyl 
operator A 4 = V 4 — 2 H 2 \7 2 in dS spacetime [9j which is a conformally covariant operator 
only when choosing M 2 = 2 H 2 . The HZ spectrum obtained by Fourier transforming the 
propagator in dS spacetime was confirmed by computing the power spectrum directly in 
momentum space. Also, the scale-invariant tensor spectrum generated during dS inflation 
could be found from the conformal gravity of ^J—gC 2 [TO] . This suggests strongly that 
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a fourth-order derivative theory related to conformal symmetry provides the HZ scale- 
invariant spectrum in whole dS spacetime. 

In this work, we wish to look for another model which may provide a HZ scale-invariant 
spectrum. This would be the nonminimal derivative coupling HB H2 m M, US] with 
fourth-order term because this coupling provides — 3£if 2 V 2 -term naturally for dS inflation. 
If a priori arbitrary parameter £ is tuned to be 2/3, one could obtain the scale-invariant 
spectrum. In this case, this model is more attractive than the LW model. Although the 
computation of the slow-roll inflation (quasi-dS expansion) is promising to compare with 
observation data, we here compute power spectrum generated during dS inflation because 
the computation is more easy and intuitive than those in quasi-dS expansion. 


2 Gravity model 


We start with the gravity model whose action is given by 


Senf — Se + Snf — / d A x\/—g 


R_ 

2 K 


A + 


1 


2 M 2 


£G' / ^0<9 i '0 - (V 2 0)‘ 


( 1 ) 


where k = 8nG = 1/Mp, Mp being the reduced Planck mass, and M 2 is a coupling 
parameter with dimension of mass squared. Here — Rxg w /2 is the Einstein 

tensor and £ is a coefficient to be adjusted as |. A is introduced as a positive cosmological 
constant to obtain dS spacetime, instead of potential V. We note that the scalar 0 has 
dimension of mass. The former term in the last parenthesis is the nonminimal derivative 
coupling (NDC) term [11, 12,[121E, 15] and the latter denotes the fourth-order kinetic (FK) 
term. We note that the former generates no more degrees of freedom (DOF) than general 
relativity canonically coupled to a scalar field, while the latter generates a new scalar DOF 
because it is a fourth-order derivative term. It may be noted that the former generates 
healthy (ghost-free) higher derivative terms, while the latter is unhealthy (ghost) higher 
derivative terms. At the first sight, the combined action Snf seems to be an unbalanced 
scalar model. For reference, we introduce the LW model given by mmm 


Slw = 


1 

2 


d A xy/^g — (V 2 0) 2 


( 2 ) 


where the first term denotes the canonically coupled (CC) term. 

Now we wish to explain the cosmological relevance of our model Snf- The NDC was 
introduced firstly by noting that the friction is enhanced gravitationally at higher energies 


3 







in the slow-roll inflation [14] • Actually, the NDC makes any potential adequate for inflation 
without introducing dangerous higher time derivative (ghost state) [13]. The NDC flattens 
a steep (non-flat) potential effectively as well as it increases friction. Also, we note that this 
kind of coupling could be view as either a subgroup of Horndeski theory [T9] or the most 
general scalar-tensor theory which is regarded as a ghost-free theory [20J. These theories 
were used to construct theoretically consistent model of inflation [2T]. Here, we wish to 
explore the role of NDC in the dS inflation by reminding that G^ takes the simple form 
of during dS phase. 

At this stage, we would like to mention the inclusion of FK term. In this case, one 
immediately ask the question of “can the model (JTI) be a consistent theory?” because the 
FK term generates the ghost in Minkowski spacetime. The existence of the ghost may imply 
the inconsistency of the model. For simplicity, we consider the LW theory (j2]) because it 
includes the FK term too. It is worth noting that this model was employed successfully 
to derive a bouncing universe to avoid singularity and a scale-invariant spectrum [ IB] . 
This was possible to occur because (J2]) can be rewritten by two fields: one is a normal 
scalar and the other is a LW (ghost) scalar <fi. The LW perturbation theory is a higher 
derivative theory and thus, it contains propagators with wrong sign residues about the 
new poles. Lee and Wick have provided a prescription for handling this issue PH HZl- 
The LW-particles associated with these new poles are not in the spectrum, but 0 decays 
to ordinary degrees of freedom instead. Their resummed propagators do not satisfy the 
usual analyticity properties since the poles are on the physical sheet. They have proposed 
deforming integration contours in the Feynman diagrams so that there is no catastrophic 
exponential growth as time increases. This amounts to a future boundary condition and 
thus, the LW theory violates the usual causal conditions. While the Lee-Wick interpretation 
is peculiar, it seems to be consistent, at least in perturbation theory, and predictions for 
physical observables can be made order by order in perturbation theory. 

We remind the reader that for £ = 2/3, the contribution from the ghost (mcc scalar) is 
found to be essential to obtain a scale-invariant spectrum by canceling of time-dependent 
term arising from the healthy (nunc) scalar in dS spacetime. Even though the FK term 
causes the ghost state, it is necessary to have a scale-invariant scalar spectrum in the 
dS inflation because it provides a scalar degrees of freedom other than 0. In the case of 
£ 7^ 2/3, S'nf implies two scalars: one is a healthy massless scalar and the other is an 
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unhealthy massive scalar. In this case, there is no cancellation between two scalars and 
thus, the scalar power spectrum is not positive definite (unitary). 

Varying the action ([Tj) with respect to the metric tensor g lw leads to the Einstein equa¬ 
tion 


GW + C + T™), 

where two energy-momentum tensors are given by 

M 2 T* dc = ^RV y ,<t>V v (!>-2R {v p V\ p \(j>V^(f>+^G^{V(f>) 2 -R^ o VyV 0 

- V /i V p 0V i/ V /O 0 + (V^V^) V 2 </> 

+ g^(R p °S7 P 4>\7A - i(V 2 0) 2 + i(V p V»V p V c 


( 3 ) 


( 4 ) 


and 


m 2 t^ = -v^(vV)v^-v v (W)v m 0 + ^[v p (vV)v^ + ^(v 2 0) 2 


( 5 ) 


We observe that T NDC r\_/ T tK when one disregards curvature coupled derivative terms i 


m 


On the other hand, the scalar equation for the action ([!]) is given by 

V 4 0 + = 0. (6) 

In order to see how the universe evolves with ©, we consider the spatially flat FRW 
spacetime by choosing cosmic time t as 

dsp RW = —dt 2 + a(t) 2 (dx 2 + dy 2 + dz 2 ). (7) 


Replacing A by the potential V, two Friedmann equations take complicated forms 


H 2 = 


K 

3 L 


Pndc + V + pfk 
9 H 2 


H = — 


- + V + — (?>H(f) 2 + - -H 2 <j) 2 - - 

3 12 M 2 ^ M 2 \ v w 2 v 2 

kt/3 H 2 

2i\W M 2 J 


1 


H M? + 


1 

M2 


(6H<j) 2 + 200 + 


( 8 ) 

(9) 


where the overdot Q denotes the differentiation with respect to time t. We observe from 
© that even though fourth-order derivative terms are generated in T)^ DC , there is no ghost 
state in pndc which means that any dangerous higher time-derivative is not generated. This 
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is why one favors the NDC term to enhance friction effects in the slow-roll inflation pH) : T5j . 
On the contrary, T™ generates ghost states in p FK because of fourth-order derivative term. 
Hence, the inclusion of p F K is problematic to obtain a solution of the universe evolution. 
This is a basic difference between NDC and FK in the FRW universe. 

The scalar equation is given by 

'0+6 H4> + ( 12H 2 + 6H)4> 

+ (9H 3 + 15HH+ 3H)<j) +M 2 —= 0. (10) 

d<p 

Requiring the slow-roll conditions of cf) 2 < 2 H 2 and (j) < 3H<p with neglecting all higher- 
order time derivative terms from FK term, Eqs. (J8]) and (TTTTh may be written approximately 
as 


n K 

H ~ 3 V ' 

(11) 

M 2 dV 

(12) 

3H ^-3 


which seem to be oversimplified equations. This means that it is very difficult to fold 
the corresponding slow-roll equations if one includes the FK term. However, one always 
obtains the dS solution for 0=const and V = A because two Friedmann and scalar equations 
reduce to H 2 = kA/ 3 and H — 0 with a trivially satisfied scalar equation. In this case, 
all scalar-derivative terms disappear but its perturbation will keep the nature of NDC+FK 
terms. This is a reason why we will choose the simplistic background geometry of dS 
spacetime to study cosmological implication of our action (ITj). ft is interesting to note that 
the Einstein+NDC+CC terms provides the dS form a{t ) = e m with H = (3y/£/M 2 ) -1 and 
</>(*) = e-VVW* at t = —oo [22]. 

A solution of dS spacetime to Eqs.flH]) and (JH]) together with = 0 can be easily found 
when one chooses a constant scalar 

R = 4kA, (j) = const. (13) 

Here, the curvature, Ricci, and Einstein tensors can be written by 

Rfiupcr (fjj/pfjurr !1 fm!]V[>): Rfiu 3H CJpui G/iv 3H Q^ (Id) 

with Hubble constant H = a/kA/ 3. Further, the flat-slicing of dS spacetime can be realized 
by introducing conformal time rj as 

ds^s = g,j,vdx ll dx v = a(r]) 2 [—dr] 2 + d %3 d;£'dx 3 }, (15) 
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where a{rf) is conformal scale factor expressed by 


a(p) = 


1 

Hrj' 


( 16 ) 


During dS inflation , the scale factor a goes from small to a very large value like af/a,i ~ 10 30 
which implies that the conformal time 77 = —1 /aH{z = — kr /) runs from — 00 ( 00 ) [the infinite 
past] to 0~(0) [the infinite future]. However, the termination of inflation is not controlled by 
the cosmological constant and thus, the inflation continues eternally. In other words, there 
is no mechanism to exit from inflation. Even though the dS inflation has such a handicap 
in compared to slow-roll inflation (quasi-dS inflation), we choose the dS inflation because 
our model (JT]) seems unlikely to provide the slow-roll inflation. 

The dS SO(l,4)-invariant distance between two spacetime points and x ,fl is defined 
by 


Z(x,x') = 1 


iv~v') 


A2 


J 12 


47777' 


= 1 


[X 


X 


A2 


47777' 


(17) 


since Z(x,x') has the ten symmetries which leave the metric of dS spacetime invariant. 
Here (x — x') 2 is the Lorentz-invariant flat spacetime distance. 


3 Scalar fourth-order propagator 


One begins with general perturbed metric with 10 DOF 


ds 2 = a(r]) 2 


(1 + 2 T)d? 7 2 + 2Bidridx l + (Sij + hij)dx l dx 2 


where the SO(3)-decomposition is given by 


(18) 


Bj — diB + 4y, hij — 2 < &5ij + 2dijE + diEj + djEi + h^j (19) 

with the transverse vectors d, x I> 1 = 0, diE 1 = 0, and transverse-traceless tensor dih l i = 
h — 0. To investigate the cosmological perturbation around the dS spacetime we 

might choose the Newtonian gauge as B = E = 0 and E* = 0. Under this gauge, the 
corresponding perturbed metric with 6 DOF and perturbed scalar can be written as 


ds 2 

0 


a ( 7 ?) 2 — (1 + 2T)d?7 2 + 2 x k i dridx l + <j (1 + 2Q)5ij + h tJ ^dx l dx 3 
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( 20 ) 

( 21 ) 









Now we linearize the Einstein equation (J3]) around the dS background to obtain the cos¬ 
mological perturbed equations. It is known that the tensor perturbation is decoupled from 
scalars. The tensor equation becomes a tensor form of a massless scalar equation 


SRfj, v (h) — — 0 —> — 0 . 


( 22 ) 


We mention briefly how do two scalars T and d>, and a vector dq go on. The linearized 
Einstein equation requires d' = — d? which was used to define the comoving curvature 
perturbation in the slow-roll inflation and thus, they are not physically propagating modes 
in dS spacetime. During the dS inflation, no coupling between {dq d>} and p occurs since 
0 = 0 implies ST^ C = = 0. Furthermore, the vector dq is not a propagating mode 

in the ENF theory because it has no kinetic term. Hence, we have the tensor hij with 2 
DOF propagating in dS spacetime. 

It would be better to find the scalar power spectrum by making Fourier transform of 
propagator in dS spacetime. First of all, we consider the £ = 2/3 case. In this case, the 
perturbed scalar equation takes the form 


vy+= o -> 5 . 2/3 vy - 27 / 2 vy = a = o, 


(23) 


where A 4 is just the Weyl operator (conformally covariant fourth-order operator) in dS 
spacetime [9] [23] . This is the main reason why we have introduced our action of S'enf ©■ 
Actually, the non-degenerate fourth-order equation can be factorized as 


V 2 (V 2 - 2 H 2 )ip = 0, 


(24) 


which implies two second-order equations for nunc and mcc scalars 


vV mmc) = 0, 

(V 2 - 2H 2 )^ mcc) = 0, 


(25) 

(26) 


where the solution to (1241) is given by 99 = + y/ mcc ). For simplicity, ( i ) denotes two 

cases: (1=1) for nunc and (1=2) for mcc. We emphasize that a choice of £ = 2/3 leads to a 
mcc scalar. Otherwise, one has a massive scalar propagating on dS spacetime. 


Expanding ip® in terms of Fourier modes (f)^ (rj) 



(27) 
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Eqs. (l25|) and (126|) become 


/ d 2 2d 

V dz 2 z dz 
/ d 2 2d 


+ 1 + 


+ 1 
2 


4 2) 


\dz 2 z dz z 

with z = —rjk. Solutions to (l28]i and (l29]i are given by 


0, 

0 


4 2) 


= ci(* + z)e 1 ' 


= c^ize , 


(28) 

(29) 


(30) 

(31) 


where c 4 and c 2 are constants to be determined. These will be used to compute the power 
spectrum directly in the next section. 

On the other hand, the linearized equation with an external source takes the form 

M 2 


A 4 y? = -M 2 J tp ip(x) = —— J v = -D(Z(x,x'))J ip (x'), 


where the propagator is given by the inverse of A 4 as 


(32) 


D(Z(x,x')) = 


2 r 


M‘ 

2H 2 


-V 2 


-V 2 + 2 H 2 


M 2 
2H 2 


[G mmc (Z(x, x')) - G mcc (Z(x, x'))] (33) 


with the dS-invariant distance Z(x,x f ) f|T7|) . Here the propagators of mmc scalar [25J and 
mcc scalar [26] in dS spacetime are given by 

H 2 


G mmc (Z(x,x')) = 


(47r) 5 


-2 ln(l - Z) + c 0 , G mcc (Z(x, x')) = 1 1 (34) 


1 - Z 


(An) 2 1 - Z 

where the former is the dS invariant renormalized two-point function (on the space of non¬ 
constant modes), while the latter is the simplest scalar two-point function on dS spacetime. 
Substituting (1341) into (l33lh the propagator takes the form 

M 2 f 


D(Z(x,x')) = - In[1 - Z(x,x')] + |) 


(35) 


which is a purely logarithm up to an additive constant Cq. 

The scalar power spectrum is defined by Fourier transforming the propagator at equal 
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time 77 = rj' as given by 

v = 1 


(2tt) 3 
1 k 3 M 2 


d 3 r Ank 3 D(Z{^x, 77 ; x', rj))e 

'i 


" ik r . r = x - x' 


(27t) 3 47T J 
1 k 3 M 2 
(27r) 3 47T 
M 2 k 2 ro ° 


d 3 r -In 


-4?7 2 - 
d 3 r ln[r 2 ]e _ik ' r 


+ 

+ 


co A e _jkr 


2 

^-(ln[477 2 ] + |)5 3 (k) 


87T 3 


dr< r sin [kr] ln[r 2 


where the last (time-dependent) term in (l38]i disappeared, thanks to 

= ™ = 0 . 

sin 0 


(36) 

(37) 

(38) 

(39) 

(40) 


We may use Cesaro-summation method to compute a logarithmically divergent integral 
m m- For this purpose, we note that the integral of J c “ f{x)dx is Cesaro summable, if 


(C, a) — hm f dx ^1 — — j f(x) 


(41) 


exists and is finite for integer a > 0. Then, (C,{3) is also Cesaro summable for any integer 

/3 > a. 

To investigate Cesaro-smnmability of the integral explicitly, we focus on f(r) = r sin [hr] ln[r 2 ] 
in (j39]h In this case, (C, a) is 


(C, a) = lim 


A —>00 


dr [l — — 


r sin[/cr] ln[r 2 


(42) 


We note that the integral (C, 0) takes the same form as (139|) . After manipulations, we have 


(C, 0) = —— lim ( kX coslfcA] — sin[/cA]) ln[A 2 / 2 ] + Si[k\] 

k 2 A->oo L v ' 


(43) 


where Si [re] denotes the sine-integral function defined by Si[x] = j^(sin[t\/t)d,t. The hrst 
two terms in (TI3|) diverge in the A —> 00 limit and thus, (C, 0 ) is not a convergent integral. 
However, the last term in (143|) is finite as it is shown in 


rk\ 


lim 

A—^00 


sinltl T 7 r 
——dt = 
t 2 


Further, (C, 1) is also not convergent since it becomes 

1 


(C, 1) = —— lim ( sin[/cA] ln[A 2 / 2 ] + SiffcAlV 

k 2 A—>-oo V / 


(44) 


(45) 
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where the first term diverges in the limit of A —» oo. On the other hand, for a = 2, the 
corresponding integral (EI2|) has a finite value 

(c ’ 2) = llb“ ( Si[a 0=-| (46) 

where we used dSD. Also, we have checked that in the limit of A —> oo, 

(C,2) = (C,3) = (C,4) = (C,5) = --- . (47) 


This implies that (C,/3) is also Cesaro summable for any integer (3 > a = 2. 
Considering (j39j) together with (1461) . the scalar spectrum takes the form 


V = 


M 2 
8tt 2 ' 


(48) 


For M 2 = 2 H 2 , it leads to an exactly scale-invariant spectrum 

Ml ) 2 


(49) 


This can be easily checked by noting that for M 2 = 2 H 2 and £ = 2/3, S'nf in ([I]) becomes 
the Lee-Wick scalar model in dS spacetime [8j - 

If one uses (fo3l) instead of (j55j) with M 2 = 2 H 2 , after Fourier transforming it, its 
spectrum is computed to be 


Vi = V n 


V mcc — 


H \2 
2tt) 


1 + ikr]) 2 — ( kr)Y 


~ (I) 2 - 


(50) 


which is a unitary scale-invariant spectrum. 

Finally, we would like to mention the unfixed £ case whose fourth-order equation is 
given by 

VV + = 0 ->■ VV - m 2 VV = = 0 (51) 

with the mass squared 

m\ = 3£H 2 . (52) 


Its propagator takes the form 
M 2 r 1 


D[Z(x, x'), m|] = 


2i/ 2 


-V 2 - 


V 2 + m 2 


•2 r 


2^2 


G mmc (^(x, x')) - G[Z(x, x'), 


(53) 
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In general, a massive scalar propagator G[Z;m|] must depend on the SO(l,4)-invariant 
distance Z(x,x') (TT7T) which has ten symmetries in dS spacetime. Also, it should satisfy 
the scalar wave equation for x ^ x' (Z ^ 1) 


i + sa-azjA 


2 

H 2 . 


G[Z(x, x'); m|] 


= 0 


(54) 


whose solution is given by the hypergeometric function 


G[Z; m|] - Y^3-r(— + ^)r(- - z^) 2 A] - + z/ ? , - - z/ 5 , 2; Z 





(55) 


At x = x', (I55j) is correctly normalized for (154)) to give <5 4 (x, x') with unit weight. We observe 
that > 0 —>■ 0 < £ < |. The negative £ is not allowed because it gives us tachyon. Here 
we note that G[Z, m 2 = 0] —* G mrac (Z) and G[Z, m| =2 / 3 = 2i/ 2 ] —» G mcc (Z). However, it 
is a formidable task to obtain its scalar power spectrum by Fourier transforming (155]) at 
rj — rf because it is difficult to Fourier transform G[Z(x, rj] x', rj), m|] with Z = 1 — 

In the next section, we could compute the scalar power spectrum for arbitrary £ < 3/4 
directly. 


4 Scalar spectrum 


In order to compute scalar power spectrum directly, we have to obtain the second-order 
bilinear action. Making use of the Ostrogradski’s formalism, one may rewrite the fourth- 
order bilinear action <5 Snf obtained by bilinearizing (CQ) as the second-order bilinear action 


xc( 2 ) _ 
uo NF — 


2 M 2 


cTx 


a 2 H 2 


— 4 )a 2 — 3£d l pd l p'j — ^(c/) 2 — 2 diCxd l ot + d 2 pd 2 p 


+4 aHaa' — 4aHad 2 p) + 2M 2 \(a — p') 


(56) 


where a = p' is a new variable to lower fourth-order derivative down and A is a Lagrange 
multiplier. In (1561) . the prime (') denotes the differentiation with respect to rj. 

From the action (|56p . the conjugate momenta are given by 

= -A, 7T a = + 2 aHa). (57) 

Varying the action (156]) with respect to A and a leads to equations 

a = p', A = ( ft// + °?H 2 ( 3£ — 2 )p' — 2d 2 a + 2aHd 2 p S J. (58) 


12 















Plugging (T58|) into f[5T]) . the conjugate momenta are given by 

= IvP (?'" + 0,2 H2 ( 3 ^ ~ 2 V - 2(9 V + 2 aHd 2 ip ), TT a = + 2aHip'). (59) 

The canonical quantization is accomplished by imposing two commutation relations 


x ), x ')] = ^( x - x '), [ & (V , x ), ^a(v, x ')] = ^( x - x ')- (60) 


The held operator (p can be expanded in Fourier modes as 


00b x ) 



d 3 k 


(ak^fa) + b ^k\v))e lii ' x + h.c. 


(61) 


where 0^ and were given by (1301) and (T3TT) . When one substitutes (IfiTT) into the operator 
of 70,, d(= <p'), and 7f Q , one obtains the corresponding expressions. Plugging these all into 
(I60p . two commutation relations take the forms 


[a k , a k /] = 5(k - k'), [ 6 k, &£,] = -5(k - k'). (62) 


It is noted that two mode operators (a k , 6 k ) are necessary to take into account of fourth- 
order theory quantum mechanically as the Pais-Uhlenbeck fourth-order oscillator has been 
shown in m- We remind the reader that the unusual commutator for (6 k , b^,) rehects that 
the FK term contains the ghost state scalar [28]. In addition, two Wronskian conditions 
are found to be 


0k } { + ((3£ - 2 )a 2 H 2 + 2k 2 ) (0*%))' - 2 aHk 2 ^(v)} 

-0k ) { (0k 2) (vi) + ((3f - 2 )a 2 H 2 + 2k 2 ) ^0* (1) (r ])) - 2aHk 2 (^ 2) ( 77 )} 


(0k } ) {(0k (1) (^)) +2a#(0 k (1) (??)) } 


- (0k } ) { (0k' 2) (v)) + 2aH (0 k (2) (r/)) } 


— c.c. = —iM 2 , 


c.c. = iM 2 , 


(63) 


which will be used to fix the coefficients of solutions 0^ and 0^ to two second-order 
equations 


(V 2 


v 2 4" 

m\)pp 


0. 

0, 


(64) 

(65) 
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( 66 ) 


where m| is given by (J52]). Explicitly, the solutions are determined to be 


(i) _ M H , ■ ( 2 ) M H 3 


m 5 VW 


('i + z)e lz , 0 ^ = 


z 2 H u (z) : 


y/2k? 

where H,^ is the first-kind Hankel function and its index is given by (155)) . 
On the other hand, the power spectrum of the scalar is defined by [6j 

( 0 |^(,,x)^ta,x')| 0 )= [ 


(67) 


Considering the Bunch-Davies vacuum state imposed by 10) = 0 and &k|0) = 0, (1671) is 
computed as 

K( k ,v) = 


where we have used (l52]h Importantly, the minus sign (—) in (168)) appears because the 
unusual commutation relation (fek^kO for ghost state was used. We expect to derive the 
power spectrum (I55|) by Fourier transforming (155|) at ij = rj'. 

For £ = 2/3= 1/2), it is shown that considering 


k 3 
2^2 
M 2 
12^C 






(i+yi-MiMM)i 




( 68 ) 

(69) 


H[% M = 

the power spectrum (JHHJ) is found to be 


z -l/ 2 e i(z-n/ 2 ) 


V = 

* 8tt 2 


M 2 r k 2 
1 + 


(aH) 2 ( aH ) 2 


M 2 

8tt 2 ' 


(70) 


(71) 


For M 2 = 2 H 2 , it leads to the HZ scale-invariant power spectrum 

2 


V 






(72) 


which is just the same form as in (149)) . We emphasize that the minus sign in (17X1) is essential 
to derive a scale-invariant spectrum from the scale-variant spectrum (1691) . 

Finally, we would like to mention that for £ ^ 2/3, its power spectrum (1591) becomes 


K(k,v) = 


Ad 2 


12£vr 2 


1 + 


/ k \ 2 7r / k 
(off) ~2\^H 




Mi 


(73) 
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which is obviously a scale-variant spectrum for dS inflation. In the superhorizon limit of 
z —* 0(k <C aif), (175j) takes the form 


V^k, rj)\ z ^ 0 = 


M 2 
12£7t 2 


(IM. 

\r(3/2) 


2ci H 


3—2uc 


which is not scale-invariant. 


(74) 


5 Discussions 

We have obtained an exactly scale-invariant spectrum of scalar perturbation generated 
during de Sitter inflation from the gravity model of the nonminimal derivative coupling 
with fourth-order term. This is the case of £ = 2/3 where there was cancellation between 
healthy massless scalar and unhealthy (ghost) scalar. For £ ^ 2/3, we have obtained the 
scale-variant spectrum (175]) where there was no cancellation between healthy massless scalar 
and unhealthy massive scalar. 

The nonminimal derivative coupling term generates a healthy second-order term of 
—2 H 2 X7 2 (p for £ = 2/3, while the fourth-order term provides an unhealthy fourth-order 
derivative term. This combination provided the linearized scalar equation of A 4 <p = 0 
expressed in term of the Weyl fourth-order operator A 4 in dS spacetime. In this sense, our 
model Snf °f NDC with FK term is more promising than the LW scalar model ([2]) where 
the NDC term is replaced by CC term of —V 2 <p. 

Now we explain a consistency of our perturbation theory in dS spacetime. The dS- 
invariant correlation function (155]) was constructed to possess full conformal invariance. We 
note that our propagator (1551) is a logarithmic correlator with the zero conformal weight. 
This means that our model (PQ) with £ = 2/3 may be considered as a minimal model 
which is allowed by unitarity of quantum field theory in dS spacetime. Our fluctuations 
which are similar to the anomaly scalar fluctuations [231121] can give rise to the HZ scale- 
invariant power spectrum (1491) in a fully SO(3,1) conformally invariant way on the dS 
horizon. This implies that that (1551) can give the unitary scale-invariant power spectrum 
(1501) in dS spacetime. However, our model with £ ^ 2/3 has led to the scale-variant 
spectrum (175]) which violates the unitarity. 

The HZ scale-invariant spectrum was obtained from Fourier transforming the fourth- 
order propagator (1551) in de Sitter spacetime. Taming a logarithmic IR divergence by making 
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use of the Cesaro summability technique, we arrived at the power spectrum of (H/2ir) 2 . 
Importantly, this HZ spectrum was also recovered by computing the power spectrum in its 
momentum space directly. In obtaining the power spectrum, we have used the Ostrograd- 
ski’s formalism and quantization scheme of Pais-Uhlcnbeck fourth-order oscillator. Hence, 
it is argued that our model (OQ) can provide a consistent perturbation theory on the dS 
background. Here, we remark that in the case of M 2 = 1, one has a scalar with zero mass 
dimension in (CQ). This corresponds to the case (7.3) in [24]. In this case, one has a constant 
power spectrum of V = 1/(87t 2 ) which seems to be trivial. 

Finally, we would like to mention the ghost issue because the FK term generates ghost 
problem. It is well known that NDC term does not contain a ghost due to the use of the 
Einstein tensor as a coupling function, while other NDC terms could contain ghosts m 
We discuss this issue by separating the dS background from the perturbation around dS 
background. This suggests strongly that one has to distinguish higher-time derivative 
terms (ghost) from ghost state (negative-norm state). As was observed from (JSJ), it is evi¬ 
dent that pndc does not contain any dangerous higher-time derivative terms (ghost), while 
Pfk involves fourth-order time derivatives (ghost). In the case of choosing dS background 
(H= const, V = A, </>=const), all time-derivatives of H and </> disappear. This implies that 
we never worry about the ghost issue on the background evolution. In the case of slow-roll 
inflation (quasi-dS), however, the ghost problem arises due to the presence of pfk- Going 
back to the perturbed scalar equation (l5Tj) . its fourth-order propagator (153|) shows the ghost 
state (negative sign in front of G[Z(x,x'),m 2 ]). However, for £ = 2/3, its propagator re¬ 
duces to the logarithmic form (j55j) in dS spacetime which shows no ghost states even though 
minus sign is present. By Fourier transforming it at equal-time leads to the scale-invariant 
spectrum (l49j) for M 2 = 2 H 2 which is free from ghost state (negative-norm state). This was 
confirmed by computing the power spectrum directly: (173]) — > (Eli. Even though the FK 
term gives rise to ghost state apparently, a special fourth-order propagator (17TT) could be 
represented by a difference of green functions between nunc (massless minimally coupled) 
scalar and mcc (massless conformally coupled) scalar propagating in dS background. The 
former corresponds to the dS-invariant renormalized two-point function, whereas the latter 
represents the simplest scalar two-point function in dS spacetimes. This is closely related to 
the fact that A 4 in (1241) becomes the Weyl operator (conformally covariant fourth-order op¬ 
erator in dS Spacetime). This difference indicates the power spectrum (l50j) which is clearly 
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ghost-free (positive-norm state) and scale-invariant. However, it seems unlikely to obtain 
the ghost-free propagator from the same operator of <9 2 (<9 2 — m 2 ) in Minkowski spacetime 
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